Stress- and temperature-dependent hysteresis in shear modulus of sohd hehum 
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The shear modulus of solid *He below 200 mK exhibits an unusual increase, the characteristics 
of which show remarkable similarities to those of the period reduction in torsional oscillator ex- 
periments. We systematically studied the drive strain and temperature dependence of the shear 
modulus at low temperatures. The hysteretic behavior depends strongly on the drive and cooling 
history, which can be associated with the thermally assisted Granato-Lucke dislocation theory. The 
phase diagram of the shear modulus is constructed on the basis of the emerging hysteretic behavior. 



An anomalous increase in the shear modulus at low 
temperatures was recently observed in solid helium 
The anomaly exhibits temperature, frequency, and '^He 
concentration dependences resembling those of the non- 
classical rotational inertia (NCRI) in torsional oscillator 
(TO) experiments 0, Si- These remarkable similarities 
have attracted theoretical and experimental attention 
0, @1- For instance, a dislocation vibration model for 
the NCRI is proposed. The temperature dependence of 
the period of a torsional oscillator containing solid he- 
lium is derived from the variation in the average pinning 
length of dislocations in solid helium 6] . Other efforts to 
provide a non-supersolid explanation for the TO response 
were also motivated bv the similarities between the two 
phenomena 0, IS S S] • The Cornell group recently re- 
ported that the TO motion was controlled by the same 
microscopic excitation that is agitated by temperature 
and stress independently, which they suggested can be 
understood within the framework of a dissipation model 
in which the elastic properties change at low tempera- 
tures [9|. Nevertheless, the resonant period reduction 
found in some earlier TO measurements, including Ref. 
, can be ascribed to the change in the elastic properties 
of solid helium filling a torsion rod [13] • Thus, the similar 
response of TO to the drive and temperature agitation 
can be understood in terms of the dislocation pinning 
mechanism. 

The shear modulus increase can be interpreted as the 
pinning of dislocations by ^He impurities according to 
the Granato-Lucke (GL) dislocation theory [llj. The 
original GL theory considers only the unpinning of dislo- 
cations by applied stress without a thermal-fluctuation- 
induced unbinding mechanism: a large stress reduces the 
number of pinning points by detaching ■^He impurities 
from dislocations, thus suppressing the shear modulus. 
Nevertheless, the effects of stress and temperature on 
the shear modulus are similar, so its temperature depen- 
dence can be explained by the variation in the dislocation 
loop length with progressive lengthening of freely vibrat- 
ing dislocation segments due to thermal evaporation of 
•^He impurities Moreover, the softening of the 

shear modulus is not a phase transition but a crossover 



of the thermally activated relaxation process from a stiff 
(pinned) state to a relaxed (unpinned) state. 

The discrepancy in the effects of the drive stress and 
temperature appears when the hysteretic behaviors dur- 
ing stress and temperature scans are considered. Strong 
hysteresis is one of the most remarkable characteristics 
found in both shear modulus and TO response measure- 
ments. When a solid sample was cooled at a high TO 
amplitude, the measured NCRI was small. The NCRI 
became large when the TO amplitude was reduced at 
low temperature. Saturation of the NCRI was then de- 
tected below the critical amplitude. However, the NCRI 
did not increase until a substantially larger TO ampli- 
tude was applied at low temperature. A similar hystere- 
sis was also reported in a number of TO experiments 
The hysteresis between the drive-up sweep and 
the drive-down sweep disappeared at temperatures above 
about 60 mK. The relaxation of the resonant period of 
TO was extended in the narrow region of the phase dia- 
gram where the hysteresis was dominant. 

Essentially the same drive- dependent hysteretic behav- 
iors appear in the shear modulus anomaly. The hysteresis 
in the shear modulus was attributed to the motion of dis- 



locations (17[ . When a solid helium sample was cooled at 
high amplitudes, the attachment of ^He impurities to dis- 
locations was hampered by their rapid motion. When the 
drive was reduced at low temperatures, the binding of the 
•^He impurities pinned the dislocations and enhanced the 
modulus. Impurity binding reduced the pinning length 
of dislocations, so larger stresses were necessary for un- 
binding the ^He impurities from dislocations. On the 
other hand, no clear thermal hysteresis appeared in the 
temperature scan. Although a similar agitation effect of 
temperature and amplitude (stress) was seen in TO (or 
shear modulus) measurements Q, the role of the drive 
stress and temperature cannot be equivalent, and the 
true nature of these agitations is not yet clear. Thus, 
it is important to investigate the thermal and mechan- 
ical hysteresis systematically to understand the subtle 
difference in the effect of drive stress and temperature 
agitations on the shear modulus. 

Here, we present systematic studies of stress- and 
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FIG. 1. (a) Shear modulus measured under constant stress 
scans at various stress values, (b) Shear modulus measured by 
sweeping stress values down and up at various temperatures. 



temperature-dependent hysteresis in the shear modulus 
of solid helium. A pair of piezoelectric shear transducers 
was employed to measure the shear modulus. The gap 
between the two piezoelectric transducers was 0.4 mm, 
as explained in a previous report 



18|. The shear mod- 



ulus is proportional to I / fV , where / is the frequency 
of the driving transducer (fixed at 1000 Hz here), / is 
the induced current at the detecting transducer, and V 
is the excitation voltage at the driving transducer. A 
solid helium sample at a pressure of 39 bar was prepared 
from commercially available ^He with a nominal isotopic 
•^He impurity of 0.3 ppm. We used the blocked capillary 
method to grow solid helium that was annealed at 1.95 
K for 3 h. 

Two different procedures were used to measure the 
stress- and temperature-dependent hysteresis of the shear 
modulus. First, the modulus was measured under vary- 
ing temperature at a constant applied strain (CS). The 
sample was cooled to 20 mK after undergoing constant 
strain at 500 mK. The shear modulus was measured dur- 
ing the warming and cooling scans. The warming scan 
data are only plotted in Fig. [Ija) since no difference be- 
tween the two CS scans was found (as also reported by 
other measurements 0); Fig. (Ha) shows the CS shear 
modulus change during the temperature scan under var- 
ious applied strains. The shear modulus exhibited the 
gradual onset of anomalous stiffening below about 250 
mK. Above this temperature, dislocations are unpinned 
from the impurities because both the thermal energy and 
strain energy overwhelm the binding energy. 

Second, the change in the modulus was measured un- 
der a varying applied strain at constant temperature 
(CT). The strain was set to 3 x 10"^ at 0.5 K, and the 
solid sample was cooled to various target temperatures. 
The CT shear modulus was measured as the strain de- 
scended stepwise to about 1 x 10"^ (CTdn scan) and 



ascended subsequently to the original value in discrete 
steps (CTup scan). 

Fig. [T{b) shows the change in the CT shear modulus 
during the drive sweep at various temperatures. The val- 
ues remained almost unchanged during cooling under a 
high strain of about 10~^. The subsequent decrease in 
the strain at low temperatures enhanced the CT modu- 
lus, which became saturated at sufBciently low strains. 
The shear modulus showed two different characteristic 
behaviors depending on the sample temperature when 
the strain was subsequently increased. The CT shear 
modulus obtained with increasing strain (CTup) did not 
trace that with decreasing strain (CTdn) below 70 mK, 
whereas the modulus was reproducible above 70 mK. The 
hysteresis was more prominent at lower temperatures, 
and the CTup shear modulus was not suppressed from 
the low-temperature saturated value under a stress of al- 
most 10 Pa. 

Fig. [2] shows the reconstructed shear moduli of various 
stresses measured during three thermodynamic trajecto- 
ries: CS, CTup, and CTdn. The shear moduli measured 
with stresses smaller than 0.2 Pa collapse onto a single 
curve. However, the three shear moduli obtained with 
stresses larger than 0.2 Pa do not agree at low tempera- 
tures, although the final stress and temperature are the 
same. 

Above the onset temperature. To, dislocations are com- 
pletely unpinned. The onset of the pinning process is ex- 
ceedingly gradual, so it is difficult to determine Tq. We 
denoted To as the temperature at which the shear mod- 
ulus deviates from the high-temperature linear depen- 
dence. To shows a monotonic drive stress dependence; a 
larger applied stress interrupts the binding of impurities 
on dislocations more violently and consequently reduces 
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FIG. 2. Shear modulus with various thermal and mechan- 
ical trajectories as a function of temperature. Color codes 
represent different protocols (thermodynamic histories). 
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FIG. 3. Shear modulus in solid '^He at various temperatures. 
Red squares are values under decreasing stress (CTdn) and 
subsequent increasing stress (CTup). Black squares were ob- 
tained during CS scan. 



the onset temperature. 

Below To, dislocations become partially pinned. The 
shear modulus in this partially pinned region is inter- 
mediate between the low-temperature saturation value 
of a strongly pinned (un-relaxed) state and the high- 
temperature value of a relaxed state. In addition, the 
shear moduli under different procedures coincide at high 
temperatures, whereas a history-dependent discrepancy 
appears at low temperatures. Two more characteristic 
temperatures related to the inherent discrepancy were 
identified. As shown in Fig. [21 the magnitude of the 
CS modulus always lies between that of the CTup and 
CTdn moduli. Above Th2, ah the thermodynamic tra- 
jectories give the same shear modulus. When the sample 
was cooled below Th2 , the CTdn modulus first separated 
from the other two moduli. The CS and CTup moduli did 
not diverge until the sample temperature was below T^i ■ 
Namely, all three shear moduli measured under various 
thermal trajectories show different values below T^i- 

The discrepancy among the shear moduli with vari- 
ous histories is also pronounced when they are plotted in 
the shear modulus-stress plane, as shown in Fig. |31 All 
the shear moduli with different histories are essentially 



identical below ah2 , which is a characteristic stress value 
measured at a constant temperature in the stress scan 
and corresponding to Th2 measured at the same stress in 
the temperature scan. Disagreement emerges when the 
applied stress increases. Between ah2 and ahi, the CTdn 
shear modulus deviates from the CTup and CS moduli, 
which agree. Finally, the CS modulus deviates from the 
CTup modulus above a^i , and all three moduli have dif- 
ferent values. The CS and CTdn shear moduli become 
identical again at a high applied stress, which eventually 
detaches impurities from dislocations and suppresses the 
CS modulus. Even under high stresses, the CTup shear 
modulus remains different from the two other moduli at 
sufficiently low temperatures. 

To better understand the effect of stress and tem- 
perature, we extracted the characteristic tempera- 
tures/stresses at given stresses/temperatures and plot- 
ted the hysteresis map in the temperature-stress plane, 
as shown in Fig. 2) We subtracted the CS modulus 
from the two CT curves and extrapolated the difference 
to zero to extract three characteristic temperatures. The 
set of characteristic points separates the phase diagram 
into two primary regions: the fully unpinned and par- 
tially pinned regions. The latter is further divided into 
hystcretic and non-hysteretic regions, where the subtle 
difference between stress and temperature effects is pro- 
nounced. Below the boundary, the distinction between 
the CT and CS shear moduli disappears, whereas the CT 
modulus differs from the CS modulus above the bound- 
ary. The fully pinned state is not available until dis- 
locations are strongly pinned by ^He impurities at low 
temperatures and low stresses. 

One may try to explain the hysteresis using the disloca- 
tion vibration model [1] . The hysteresis can be associated 
with the shortening or lengthening of vibrating segments 
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FIG. 4. Stress-temperature hysteresis map of solid '''He ob- 
tained at characteristic points as described in the text. 
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between two stable nodes of dislocations by ^He impuri- 
ties. During the CTdn scan, impurities are bound to dis- 
locations when the dislocation loop length is shorter than 
the critical length Lc{a). When the scan is reversed to 
a CTup scan, the previously shortened dislocation loops 
require greater stress to separate the impurities, which 
leads to the hysteresis. Thus, the hysteresis appears 
with variations in the dislocation loop length, which de- 
pend strongly on the thermal and mechanical trajecto- 
ries of solid helium. However, this vibration model is 
constructed on the nave assumption that the spacing be- 
tween ^He impurities is identical, which cannot explain 
the details of the hysteresis. For example, once '^He im- 
purities become unpinned from dislocations, stresses far 
below the critical value are sufficient to detach neigh- 
boring "^He impurities on dislocations. Accordingly, the 
breakaway of ^He impurities may induce a sudden de- 
crease in the shear modulus when the stress exceeds the 
critical value. One may use the distribution of the impu- 
rity pinning length as well as that of the network pinning 
length to explain the gradual decrease in the CT shear 
modulus during the drive up-scan po| . 

A significant feature of Fig. [2] is the discrepancy be- 
tween the CS and CTdn moduli at low temperatures. 
The discrepancy is not expected when the stress and tem- 
perature have equivalent effects on the shear modulus 
via the same microscopic excitation. When the stress or 
temperature decreases, the dislocations are progressively 
pinned by "^He impurities. The pinning and unpinning 
of impurities on dislocations are simply determined by 
the competition between the pinning potential and the 
thermal or mechanical energy. Thus, if the thermal and 
mechanical energy controlled the shear modulus indepen- 
dently, the dependence of the shear modulus on the tem- 
perature and the square-root of the strain would show 
similar responses without a discrepancy 9] . 

Fig. |4] provides additional evidence that the stress and 
temperature affect the pinning of dislocations in differ- 
ent ways. The boundary between the hysteretic and no- 
hysteresis regions has a positive slope, indicating that a 
higher temperature is required at higher stress to com- 
pletely eliminate the hysteresis. This positive slope can- 
not be understood as well if the effects of temperature 
and stress are equivalent. If the shear modulus is con- 
trolled by a simple combination of thermal and mechani- 
cal agitation, the boundary would have a negative slope. 
Consequently, a complicated interplay of thermal and 
mechanical influences controls the hysteresis map, and 
it is crucial to investigate the roles of stress and temper- 
ature in order to understand the hysteresis more thor- 
oughly. 

We used the thermally assisted breakaway model to 
describe better both the temperature and drive depen- 
dences of the shear modulus change in solid helium at 
2^ . Assuming that the interaction 



low temperature [21 



between an impurity atom and a dislocation is given as a 



simplified Cottrell force 23|, the potential energy of the 
system consists of three terms: a dislocation tensile line 
energy and an impurity-dislocation binding energy that 
likely lead to a pinned state, and the strain energy, which 
favors an unpinned state [21[ . The total potential energy 
of the system has a double-well shape with two minima 
as a function of the distance between the impurity and 
the dislocation line: a pinned state at a short distance 
and an unpinned state at a longer distance. Thermal en- 
ergy can assist the transition between the two minimum 
states. The shape of the potential energy is remarkably 
susceptible to the applied stress At low stresses, 

only one energy minimum appears as a strongly pinned 
state, and at high stresses, only the unpinned state is 
energetically stable, indicating a fully unpinned state. In 
the intermediate range of stresses, two energy minima 
appear; both pinned and unpinned states are possible. 

This intriguing model very clearly explains the hys- 
teresis of the CT shear modulus shown in Fig [1] If solid 
helium is cooled under a high applied strain, the shear 
modulus remains small because the unpinned state is sta- 
ble under a high stress. As the stress decreases, the sec- 
ond minimum (unpinned state) diminishes, and the sys- 
tem shifts progressively to the pinned state. When the 
stress is then increased at low temperatures, the second 
minimum is restored. However, the system cannot tran- 
sition to an unpinned state if the thermal energy is insuf- 
ficient to overcome the potential barrier between the two 
states. This irreversibility produces the discrepancies be- 
tween the CTdn and CTup moduli at low temperatures. 
Namely, the hysteresis in the CT modulus can be associ- 
ated with the stress-dependent effective potential barrier 
between the two minima. At high temperatures, the ab- 
sence of the hysteresis can be explained by the presence 
of sufficient thermal energy, which allows the system to 
transfer to an unpinned state. 

This framework can explain why stress and temper- 
ature have different effects. The applied stress changes 
the shape of the interaction potential; thus, the preferred 
direction of the process reverses with the sweeping direc- 
tion. In contrast, the thermal effect enhances the transi- 
tion probability between the two states without a direc- 
tional preference. The positive slopes of T^i and Th2 as a 
function of applied stress in the phase diagram reflect the 
different roles of temperature and stress. A high applied 
stress enhances the imbalance between the two states in 
the double-well potential, which causes the difference in 
the transition probability: the unpinning transition is fa- 
vored over the repinning transition because of the raised 
potential barrier from an unpinned to a pinned state. 
In this partially pinned state, the initial population in 
the two states and the thermomechanical trajectory that 
modify the population strongly influence the flnal state. 
Thermal fluctuation is the essential ingredient that neu- 
tralizes the imbalance. Therefore, higher thermal energy 
is necessary for compensating the raised potential barrier 
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in the repinning transition. 

The subtle difference between the roles of stress and 
temperature is more pronounced when the CS modulus is 
compared with CT modulus. Above Th2^ the sufficiently 
larger thermal energy compared to the transition barriers 
allows an unobstructed transition between the two states. 
Thus, the shear modulus is independent of the stress and 
temperature history. However, the population in the pin- 
ning region during a CTup or CS scan is greater than 
that during the CTdn scan between Thi and Th2 because 
most dislocations remain in the unpinned state during 
the CTdn trajectory. The increase in the drive stress 
deepens the unpinned state in the double-well potential, 
so the repinning transition from the unpinned state to 
the pinned state is essentially restricted below Th2 owing 
to the lack of thermal energy. Moreover, the popula- 
tion cannot be augmented simply by lowering the stress 
in certain temperature ranges. This restraint causes the 
strongly history-dependent behavior of the modulus. 

In contrast, if the passage is chosen to follow the CTup 
procedure, most dislocations stay in the pinned states. 
The lack of a thermally assisted crossover significantly 
suppresses the transition probability to the unpinned 
state. This observation explains why the CTup modu- 
lus always exceeds the CS modulus below Thi- 

The CS shear modulus manifests the change in the 
populations of the two states as the thermal fluctuation 
decreases and is intermediate between the two CT mod- 
uli. Above To, the modulus remains at the value of the 
unpinned state until the systems thermal fluctuation is 
reduced sufficiently below the typical binding energy of 
•^He impurities and dislocations. 

The shape of the shear modulus phase diagram is sim- 
ilar to that of NCRI It is probably reasonable to 
say that the phenomena are connected, considering their 
similarities in many properties, including the phase di- 
agram resemblance. One plausible explanation is that 
the fluctuating mass currents caused by mobile disloca- 
tions induce phase fluctuations and destroy the quan- 
tum coherence for supersolidity. Thus, the pinning of 
dislocation lines may be necessary for the appearance of 
NCRI. All the peculiar properties in the low-temperature 
anomaly of solid helium also stem from this association. 
However, the stress values in the NCRI phase diagram 
are much smaller than those of the shear modulus [l8| . 
Furthermore, an experiment in which DC rotation was 
superposed on the shear modulus and TO measurements 
demonstrated a clear contrast between the two phenom- 
ena, suggesting that they cannot have the same micro- 
scopic origin [24)| . Therefore, the underlying relationship 



between the shear modulus and the TO anomaly is very 
intricate, and a direct probe of both anomalies, such as 
simultaneous measurements of both phenomena, is cru- 
cial. 
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